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In a two level atom, real-time quantum state holography is performed through 
interferences between quantum states created by a reference pulse and a chirped 
pulse resulting in coherent transients. A sequence of several measurements allows 
one to measure the real and imaginary parts of the excited state wave function. 
These measurements are performed during the interaction with the ultrashort laser 
pulse. The extreme sensitivity of this method to the pulse shape provides a tool for 
electric field measurement. 

I. INTRODUCTION 

The effect of laser pulse shape on a quantum system is related to the nature of the inter- 
action. For a linear response of the system (one-photon transition in the weak field regime), 
the final outcome depends only on the spectral component at the resonance frequency and 
is therefore independent of the pulse shape, and particularly of the spectral phase [lj]. This 
explains for instance why signals equivalent to wave-packet interferences could be observed 
with incoherent light as well as with ultrashort pulses However, the temporal evolution 
towards the final state may depend strongly on the pulse shape. A straightforward illustra- 
tion of this statement is the non-resonant interaction which leads to transient excitation of 
the system, but to no final excitation. In the absence of predesigned control mechanisms 
only a dosed ,oop Scll e m e m ay be e mp l„yed to Had efficient pulse shapes = 
The outcome of many different shapes is fed back into an algorithm that iteratively opti- 
mizes the excitation shape without insight into the physical mechanism that is triggered by 
a particular shape. 

In contrast the effect of shapes on small systems can be systematically studied within an 

, This open-loop approach is well adapted to these systems for 
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which theoretical predictions are reliable. It consists of reaching a specific goal (manipulation 



of the temporal response of a system excited by a light pulse) without any experimental feed- 
back. Physical analysis of the process allows one to predetermine the theoretical pulse shape 
which leads to the desired result. It is then implemented experimentally. 

In this article, we describe manipulation of Coherent Transients (CT) in an open loop 
approach. These CT are oscillations in the excited state population resulting from the inter- 
action between a two-level system and a weak chirped pulse. The shape of these oscillations 
is extremely sensitive to slight changes in the pulse shape 3, Q|- Two previous letters 
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have shown that their high sensitivity provides a new scheme for quantum state 
measurement and electric field reconstruction. This article presents in details the works and 
calculations corresponding to these results. First we recall the coherent transients and how 
to manipulate them. Then the quantum state measurement reconstruction is presented in 
detail. In particular several schemes are discussed. Then the experimental set-up and several 
previous feasibility test are described. Finally the results are presented and discussed. 



II. COHERENT TRANSIENTS PRINCIPLE 



The CT result from the interaction of a two-level system (\g) and |e}) with a chirped 
pulse E pu {t) of carrier angular frequency ui pu close to resonance (u) pu ~ u) eg ). The transient 
excited state population is probed towards the |/) level in real time by a second ultrashort 
pulse E pr {t) which is Fourier transform limited and short compared to the characteristic 
features of E pu {t). Its frequency is close to resonance (ujf e ). The fluorescence arising from 
the \ f) state is then recorded as a function of the pump-probe delay r (cf Fig.HJ. The probe 
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Fig. 1 - Excitation scheme. 



pulse provides access to the temporal evolution of the population in |e), produced by the 
pump beam. The result of the interaction is described by first order perturbation theory, 
and the fluorescence is proportional to 
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in the case of a simply chirped pulse E pu (t) = e~ (t/T P u)2 e~ i(Wpil ' +apil ' 2) . Here 5u eg 
is the resonance mismatch, r pu the pulse duration and a pu the chirp rate. A quadratic phase 
appears in the integral giving q e (t) (Eq. EJ, leading to oscillations of the probability |a/(r)| 2 
as already demonstrated (cf Fig.EJ) • These strong oscillations result from interferences 

between the probability amplitude excited at resonance and after resonance. 




2 4 6 

Delay (ps) 



Fig. 2 - Experimental Coherent Transients on Rb (5s-5pi/2 at A = 795 nm), for a chirp of 
—8. 10 5 fs 2 (dots) and the corresponding simulation obtained by numerical resolution of the 
Schrodinger equation (solid line) [ifj. Inset : Theoretical excited state amplitude drawn in 

the complex plane. 



The CT phenomenon is better understood by examining the behavior of a e (t) in the 
complex plane as displayed in the inset of Fig. [21 The probability amplitude follows a Cornu 



spiral starting from the origin. Three regions can be distinguished. The two spirals result 
from contributions before (I) and after (III) resonance for which the quadratic phase varies 
rapidly. The intermediate region (II) corresponds to the passage through resonance where the 
phase is stationary. It provides the main contribution to the population. The two spirals, 
although similar, result in totally different behaviors of the population. The first one (I) 
winds round the origin with an increasing radius. The resulting probability increases thus 
slowly and regularly and remains small. After resonance (III), a second spiral winds round 
the asymptotic value leading to strong oscillations of the population. 

We show in the next section how a modification of the excitation scheme provides the 
possibility to observe oscillations due to the first part of the pulse. 



III. QUANTUM STATE MEASUREMENTS 



A. principle 

The CTs are extremely sensitive to tiny phase modifications of the pump pulse Q|- 
Therefore, they can provide detailed information on the exciting pulse and simultaneously 
on the excited quantum state. However, although sensitive to phase effects these CTs give 
access to the excited state probability |a e (r)| 2 whereas the probability amplitude is necessary 
to achieve a complete measurement of the electric field. Moreover, the oscillations are only 
produced by the second part of the pulse (after resonance) To overcome these limitations, 
we propose a new excitation scheme based on a two pulse sequence with a well defined phase 
relationship. The pump pulse is written as 

E pu (i) = E pui (t) + e t6 E pU2 (t) (3) 

where E pui (t) and E pU2 (t) are two replica of the initial pulse with additional spectral phase. 
These can be obtained with either a Michelson-type interferometer or a pulse shaper. The 
first pulse E pui (t) creates an initial population in the excited state. The second pulse E pU2 (t) 
is strongly chirped and sufficiently delayed so that it does not overlap with the first pulse. 
This second pulse creates a probability amplitude in the excited state which interferes with 
the initial probability amplitude created by the first pulse. 

It should be noted that the details of the shape of the first pulse are not critical. Only 
the final state reached at the end of the first pulse is involved in the temporal evolution of 



the system during the second pulse. 
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Fig. 3 - (a) Theoretical CTs scans for a geometric reconstruction : 6 = (black), 6 
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2.10 5 fs 2 . (b) Corresponding probability amplitudes (same 
color code). 



Looking at the evolution of the quantum state in the complex plane (Fig. 02(b)), one sees 
that the effect of the first pulse is to shift the starting point of the spiral so that oscillations 
due to CTs occur on the whole duration of the second pulse. Assuming a sufficient time 
interval between the two pulses to avoid any overlap, the probability amplitude induced by 
the first pulse a el (t) has reached its asymptotic value a el (oo) when the interaction with the 
second pulse starts. For a probe pulse significantly shorter than the details one wants to 
retrieve on the excited state population, the recorded fluorescence S e (r) is directly propor- 
tional to the excited state population. During (or after) the second pulse, it can be written 
as 



a e i(OOl 



e l6 a e2 (r)\ 2 



;1 (oo)| 2 + |a e2 (r)| 2 + 2Re [e i6 a e l(oo)a e2 \ 
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|a el (oo)| can be deduced from a measurement of S e (r) in the interval between the two 
exciting pulses. In order to determine the complex number associated to the probability 
amplitude, at least a second measurement is necessary as described in the next subsection. 



B. Reconstruction techniques 



The probability amplitude produced by the second pulse E pU2 (t) is retrieved by combining 
the CTs scans S e (r) (see Eq. EJ) obtained for different values of the programmable phase 



9. The goal here is to extract the cross term a e \{oo)a e2 {r) from a set of scans. The factor 
a e l(oo) is deduced - except for its phase - from the fluorescence observed at the end of 
the first pulse. We will show here different possible reconstruction schemes. As an example, 
we simulate the CTs corresponding to the following two-pulse sequence : the first pulse is 
100 fs long, the second one is chirped to 10 ps (2.10 5 fs 2 quadratic phase) and delayed by 
6 or 10 ps. Both pulses are resonant oo pu = u eg and have the same energy. Determining the 
real and imaginary part of the probability amplitude requires at least two equations, which 
means two CT scans with different values of 9. In this case, a system of two second order 
equations is obtained. A geometric method is used to solve it. With a third measurement, 
the quadratic term in Eq. 0]can be removed in order to obtain a linear system. 

With a set of three scans, the angles 9k(k = 0, 2) must be chosen so that the corresponding 
matrix is not singular. This is the case for instance with 9k = 2kn/3. The corresponding 
CTs are plotted in Fig. EJ From these measurements, we calculate S^r) defined as : 
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Fig. 4 - Theoretical CTs for a three-scans reconstruction : pulse sequence with (jr P u = 6 ps, 
4« = -2.10 5 fs 2 and 9 k = 2kn/3. 9 = : light grey line ; 9 X = 2n/3 : black line ; 

6*2 = 4:71 / 3 : grey line. 

S [3] (r) = \s°(r) - ii±!^W(T) - !^S^(r) = a e \{oo)a e2 {r) (5) 
3 b o 

Alternatively, the probability amplitude can be retrieved from a set of two CT measure- 
ments provided that a system of two nonlinear equations is solved. For two different values 



of 9, for example and a^O, n, we thus have to solve the two-equation system : 

S°(r) = |a el (oo)| 2 + \a e2 (r)\ 2 + 2Re [a e *(oo)a e2 (r)] (6a) 
S a (r) = |a el (oo)| 2 + |a e2 (r)| 2 + 2Re [e ia a e l(oo)a e2 (r)] (6b) 

If the second pulse is much weaker than the first one, the quadratic term in |a e2 (T)| can be 
neglected to obtain a simple linear equation system. In this case one easily obtains a unique 
solution and a = n/2 is the simplest choice. 

Generally, the non-linear equation system gives two different solutions and only one is 
physically acceptable. To easily identify this solution and separate it from the other one, we 
have developed a geometric reconstruction which is described in detail in the appendix. 



IV. EXPERIMENT 
A. Experimental set-up 
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FlG. 5 - Experimental set-up. NOPA : non colinear optical parametric amplifier; CPA 

chirped pulse amplifier 



The experimental set-up is displayed in Fig. El The 5s - 5p (P1/2) transition (at 795 
nm) is resonantly excited with a pulse sequence. The transient excited state population is 
probed "in real time" on the (5p - (8s, 6d)) transitions with an ultrashort pulse (at 607 
nm). The laser system is based on a conventional Ti : Sapphire laser with chirped pulse 
amplification (Spitfire Spectra Physics) which supplies 1 mJ -130 fs -795 nm pulses. Half 
of the beam is used for the pump pulse. The remaining seeds a home made Non-collinear 
Optical Parametric Amplifier (NOPA) compressed using double pass silica prisms, which 
delivers pulses of a few microJoule, 30 fs -FWHM pulse intensity, centered around 607 nm. 
The pump pulse is shaped with a programmable pulse-shaping device producing the pulse 



sequence, recombined with the probe pulse and sent into a sealed rubidium cell. The pump- 
probe signal is detected by monitoring the fluorescence at 420 nm due to the radiative cascade 
(8s, 6d) - 6p - 5s collected by a photomultiplier tube as a function of the pump-probe delay. 
In parallel, a cross-correlation of the pump pulse sequence is recorded. The pulse shaping 
device is a 4f set-up composed of one pair each of reflective gratings and cylindrical mirrors. 
Its active elements -two 640 pixels liquid crystal masks- are installed in the common focal 
plane of both mirrors. This provides high resolution pulse shaping in phase and amplitude 
|l7l |. This is used to generate the shaped pump pulse sequence by applying the function 




The laser is centered at resonance {ui pu = uJ e g)- 



B. Interferometric stability 

The relative stability of the two pulse sequence is a crucial point in the present experiment. 
Both the relative phase and delay between the two pulses of the pump sequence should be 
kept stable as compared to 2tt or to the optical period T Q . Experiments of wave packet 



interferences have been performed with a Michelson interferometer usee 



pulse pair. The delay was either actively [1$, |1JJ, 120] or passively 



to produce the 
2511 stabilized 
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using different techniques. The best achieved stability is better than T o /100 with a Michelson 
placed under vacuum 

Alternatively, in experiments where only the amplitude of the interference pattern is 
needed, different strategies have been developed. These are based either on periodic modu- 
lation of the delay followed by a lock-in amplifier (3, or random fluctuations applied to 
the delay followed by measurement of the resulting noise 0, S] • 

In our approach, the required stability and control of the phase and delay are naturally 
provided by the phase and amplitude pulse shaper [17]. As an illustration, we have performed 
demonstration experiments with a pump sequence consisting of two identical Fourier trans- 
form limited pulses, delayed by 3 ps. In a first example (see Fig. I0J) , the relative phase (at the 
resonance frequency) is set to and 7r for two scans of the pump - probe delay. Two cross- 
correlations measurements (Fig. UK and b) illustrate the relative position of the pulses. The 
phase shift of n does not affect these cross-correlations. In the pump-probe scan, the three 



positions of the probe pulse with respect to the pump sequence lead to : (i) No fluorescence 
signal when the probe is before the pump pulses, (ii) a constant signal independent of the 
relative phase for the probe before the pump pulses, (iii) constructive (9 = 0) or destruc- 
tive (9 = 7r) interferences for a probe pulse after the pump sequence. In the constructive 
interference case, the fluorescence signal is 4 times the signal obtained with a single pulse, 
as expected from usual interferences. 
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Fig. 6 - Experiments with a pump sequence of two Fourier Limited pulses centered at 795 
nm and a 25 fs probe pulse at 607 nm, as a function of the probe pulse delay. Pump-probe 
cross-correlations for a relative phase of (a) or tt (b). (c) Fluorescence from the 8s-6d 
states for the two relative phases (black line : 9 = ; Gray line : 9 = tt) . 



In a second experiment, the pump-control delay is set to a constant value of 267 ps and 
the relative phase is scanned (Fig.EJ). These two results illustrate both the excellent stability 
of the set-up and the control over the programmable phase. 
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Fig. 7 - Same scheme as in Fig. El but with a fixed pump - probe delay and a variable 
relative phase. Experiment (Squares) and sine fit (solid grey line). 

V. RESULTS AND DISCUSSION 

We present a series of results obtained with a sequence of two pulses generated by the 
high resolution phase and amplitude pulse shaper : The first one is close to Fourier limited 
(a replica of the input pulse) and the second one is time delayed (</4« = 6 ps) and strongly 
chirped (<fipu = — 2.10 5 fs 2 ). Their amplitude are set equal. With phase and amplitude shap- 
ing, an extra relative phase 9 can easily be added to the second pulse. The various records 
correspond to different values of 9 and are used to illustrate the two reconstruction methods 
described in Section IHIBI 

In the first example, three scans with phases separated by 2ir/3 are used : 9 = 9 + 
2kii/3 (k = 0,1,2). The quality of the reconstruction does not depend on 9 and here we 
have #0 ~ 0.7 (the reconstruction efficiency does not depend on 6*0). The scans are displayed 
in Fig. [Hl-a. As a difference to the case of a single chirped pulse (Fig. [2]) [3], the three 
regimes are now clearly visible. The oscillations are observed before resonance as well as 
after resonance. The behavior during the passage through resonance depends directly on the 
relative phase 9. A rapid increase, slow increase or slow decrease is observed resulting from 
constructive, partially constructive or destructive interferences. As expected, and similarly 
to the case of two FT limited pulses (see Fig. and [Z|) , the asymptotic value depends also 
strongly on 9. The linear reconstruction method is used. The good stabilities of the laser 
and experimental set-up allow us using directly the raw data without any adjustment. The 



excited state probability amplitude produced by the second pulse (with 9 = 9 ) is extracted 
from the three measurements and displayed in Fig|Bl-b). One observes clearly the expected 
Cornu spiral. 



Fig. 8 - a) Experimental Coherent Transients resulting from the excitation of the atom by 
a FT limited pulse (at time r = 0) followed by a chirped pulse ((/)$} = 6 ps, 
-2.10 5 fs 2 ), for three different relative phases #0 ~ 0.7 (light grey line), 9 = 9 + 2-7r/3 
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(black line) and 9 = 9q + An/ 3 (dark grey line) between the two pulses, b) Probability 
amplitude a ei 2(r) reconstructed from the three measurements presented in a), using a 
linear reconstruction and displayed in the complex plane. The Cornu spiral appears clearly. 

In the second example displayed in Fig. two scans with phases separated by n/2 : 
9 = 9 and 9 + n/2 are used for the nonlinear reconstruction. Here 9q ~ —0.8. The non- 
linear method requires determining separately the population |a el (oo)| 2 created by the first 
pulse. The plateau immediately after the end of the first pulse is used for this purpose. The 
excited state probability amplitude produced by the second pulse and extracted from the 
two measurements is displayed in Fig. E3b. The reconstructed probability amplitude is also 
displayed in Fig. [TQl in a 3D plot (real and imaginary part of the probability amplitude as a 
function of time). The projections on the various 2D planes are also displayed. The expected 
Cornu spiral is clearly seen in the complex plane. 

In previous experiments Q|- only the excited state probability was measured. Here, 
the initial preparation of a coherent superposition of |e) and \g) by the first pulse allowed 
measuring the probability amplitude in real time during its evolution in interaction with the 



FlG. 9 - a) Experimental Coherent Transients resulting from the excitation of the atom by 
a FT limited pulse (at time r = 0) followed by a chirped pulse (0p« = 6 ps, 
<Ppu = — 2.10 5 fs 2 ), for two different relative phases 6$ ~ —0.8, 6 = 6 Q + it/2 between the 
two pulses, b) Probability amplitude a e ^{j) reconstructed from the two measurements 
presented in a) and displayed in the complex plane. The Cornu spiral appears clearly. 




Fig. 10 - 3D spiral representing the time evolution of the excited state probability 
amplitude (same data as in Fig. EJ). The vertical axis represents the time. 

laser pulse. 

The two methods provide similar quality of reconstruction. The linear approach requires 
three measurements. It increases the recording time by a factor of 1.5 as compared to the 
nonlinear method. Conversely it is more robust and can be used in a wider variety of situ- 
ations (relative phase, intensity ...). A larger number of recordings could be combined in a 
linear square fit approach to improve the accuracy of the measurement. This would of course 
be at the extent of the recording time. 



Several examples of quantum phase measurements of states created by ultrashort pulses 
are based on interferences between an unknown wave function and a "reference" wave func- 
tion. These wave functions are created by a sequence of two ultrashort pulses (an unknown 
pulse and a reference pulse). The quantum state created by the unknown pulse is deduced 
either by time- and frequency- integrated fluorescence measured as a function of the delay 
(si | . or by measuring the population of each eigenstate for different values of the relative 
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phases |22j. Alternatively, the amplitude of noise resulting from random fluctuations of the 
delay is measured 0, [j3|. In another approach, the dispersed fluorescence emitted by an 
oscillating nuclear wave packet in a diatomic molecule was recorded as a function of time 
. In this case, the fluorescence wavelength - position relationship is derived from the 
Franck-Condon principle. 

In all these examples involving several excited states, either a particular selectivity is 
used to detect independently each excited state, or the delay is used to obtain a set of 
measurements which are then inverted to obtain the amplitude of each quantum state. In 
our study, only one single excited state is involved and the measurements are performed at 
the same delay. This ensures determining the temporal evolution of the quantum state. 

Our quantum state measurement method can be extended to the case of p excited states 
(\i))i=i, P of different energies ftuv Their probability amplitudes can be retrieved from a set of 
2p+l measurements in a linear reconstruction scheme. As an example, the first measurement 
can be performed with the second pulse E pU2 (t) only, providing 

S 9 (r) = J> 2 (r)| 2 (8) 

i=l 

This allows thus to remove the nonlinear contributions from the subsequent measurements. 
The remaining 2p measurements are performed with the two pulse sequence, each with a 
set of p phases (Oi,k) i= i p for k — l,2p applied at the frequencies a>j. They provide with the 
quantity 

p p p 

3 h (r) = KiMI 2 + Yl \ a ^\ 2 + 2 J2 Re [e^ fc <(°o)a i2 (r)] (9) 

i=l i=l i=l 

P 2 

As in the case with a single excited pulse, l a ii(°°)l can t> e deduced from a measurement 

i=l 

of S k (r) in the interval between the two exciting pulses. Since the phases can be chosen 
independently, it is always possible to find a set of phases for which the system of 2p lin- 



ear equations can be inverted. This would not be the case if the phases were not applied 
independently but through an extra delay r' k (giving 6^ = u} ig r'^). 

VI. CONCLUSION 

We have presented a new method to determine the real time evolution of an excited 
quantum state in interaction with an ultrashort laser pulse. 

By simple derivation of the excited state probability amplitude, it is possible to retrieve 
the electric field (phase and amplitude) of the second pump pulse E pU2 (t), provided that 
the probe pulse is well known or short enough. The possibilities offered by this technique 
are discussed in detail elsewhere It can also be used for a differential measurement to 
analyze the changes induced by inserting a material. In this last case, the requirements on 
the properties of the probe pulse are less severe. 

We thank Chris Meier for fruitful discussions. 



VII. APPENDIX 



We detail here the geometrical reconstruction used to solve the set of second-order non- 
linear equations. This latter interprets the equation system ((HJ) in terms of circles intersection 
in the complex plane. Figured (a) shows the two CTs scans for 6 = (black) and 6 = a (gray) 
used for the reconstruction (in the simulations a = vr/3). The corresponding probability 
amplitudes are plotted in Fig. El (b). In both cases, the contribution of the first Fourier- 
limited pulse is a straight line and the contribution of the second pulse is a Cornu spiral. 
The phase 9 only changes the relative orientation of the line and the spiral. At any time r, 
the CTs values S°(t) and S a (r) respectively correspond to Tq(t) and r^,(r), where r (r) and 
r Q (r) are the distances in the complex plane between the origin and the current positions on 
both spirals (see Fig. E](b)). Retrieving the probability amplitude produced by the second 
pulse corresponds to geometrically reconstructing the black Cornu Spiral in Fig. 01(b), using 
the two time dependant distances r (r) and r a {r). 

To achieve this, we mentally rotate the gray path by an angle of —a, around the starting 
point of the Cornu Spiral (1,0). We then choose this point as the new origin for the complex 
plane. These transformations preserve both angles and distances and therefore do not change 



our equation system. Figure [TT] shows the two paths after the transformations. We call 
P and P a the starting points of each path whose coordinates are known : (—1,0) and 
(— cos(a), sin(a)) respectively. Thanks to these transformations, the two Cornu Spirals are 
now superimposed and correspond to the amplitude probability we want to retrieve. The 
two distances r (r) and r a (r) can now be seen as the distances between the point a e2 (r) on 
the Cornu spiral and two reference points Pq an d P a - To geometrically reconstruct a e2 (T) 




Fig. 11 - New geometric interpretation : r (r) and r a {r) are the distances between a e2 (r) 

and two reference points Pq et P a . 

one just needs to find, for every time r, the intersection of the circle Cq (centered on P 
with a radius r (r)) and the circle C a (centered on P a with a radius r a (r)), as depicted in 
figure H21 We get two different solutions, a+ 2 (r) and a~ 2 (r). The physically acceptable one 
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Fig. 12 - Geometric reconstruction of a e2 (r). Two solutions a+ 2 (r) and a e2 (r) are 
available; The physical one starts in (0,0). 



starts in (0,0). To avoid degeneracy, the Cornu spiral should not cross the (P ,P a ) line. 



Two ways of pushing the spiral away from (P , P a ) are increasing the intensity of the first 
pulse, and reducing the angle a. However, a too small angle leads to near-degeneracy of the 
circles, increasing thus the uncertainties in determining their crossing points. Usually, the 
reconstruction works well with a first pulse at least as intense as the second one and an 
angle a in the interval [n / 4, n/2]. 
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